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Abstract. 

We present the first results of an observational programme undertaken to map the fine structure line emission of 
singly ionized carbon ([CII] 157.7409 /im) over extended regions using a Fabry Perot spectrometer newly installed 
at the focal plane of a 100 cm balloon-borne far-infrared telescope. This new combination of instruments has a 
velocity resolution of ~ 200 km/s and an angular resolution of 1 '5. During the first flight, an area of 30'xl5' 
in Orion A was mapped. These observations extend over a larger area than previous observations, the map is 
fully sampled and the spectral scanning method used enables reliable estimation of the continuum emission at 
frequencies adjacent to the [CII] line. The total [CII] line luminosity, calculated by considering up to 20% of the 
maximum line intensity is 0.04% of the luminosity of the far-infrared continuum. We have compared the [CII] 
intensity distribution with the velocity-integrated intensity distributions of 13 CO(1-0), CI(l-O) and CO(3-2) from 
the literature. Comparison of the [CII], [CI] and the radio continuum intensity distributions indicates that the 
largescale [CII] emission originates mainly from the neutral gas, except at the position of M43, where no [CI] 
emission corresponding to the [CII] emission is seen. Substantial part of the [CII] emission from here originates 
from the ionized gas. 

The observed line intensities and ratios have been analyzed using the PDR models bv lKaufman et all (I1999T) to 
derive the incident UV flux and volume density at a few selected positions. The models reproduce the observations 
reasonably well at most positions excepting the [CII] peak (which coincides with the position of 1 Ori C). Possible 
reason for the failure could be the simplifying assumption of a homogeneous plane parallel slab in place of a more 
complicated geometry. 

Key words, infrared: ISM - ISM: lines and bands - ISM: individual (Orion A) 

1. Introduction tensities of C + , C° and O lines but fail to explain the ex- 
tended nature of the [CII] emission. However a more rcalis- 

In recent years it has been established that molecular tic model suggests that the molecular clouds are clumpy in 

clouds show definite evidence of internal structures in the naturc which allows the UV radiation to penetrate deeper 

forms of clumps and filaments. The UV radiation, from and hcnce the observed emission consists of the contribu- 

external or embedded OB stars, incident on the molecular tion from y photodisso c iated clump surfaces along the 

clumps dissociates molecules and ionizes atoms within a line of sight feoster et all 11994) . The observed [CII] as 

shell of A v a 5 m on the surface of the clumps and forms wcll as the co li ne intensities scale with the density of the 

the Photon Dominated Regions (PDR). The externally PDRs and thc intensity of the UV radiation. Most of the 

heated PDR gas cools via characteristic FIR and sub-mm high spatial reso lution observational studies of [CII] line 

lines of various atomic and molecular species. The fine emission have concentrated on regions of high density and 

structure lines of [CII] 158 /xm and [OI] 63 fim typically ac- high uv radi ation intensity . Howe ver large scale surveys 

count for 1% of the cooling of the PDRs. One dimensional with C0BE JWright oi a 1. 1. Il99ll>. Balloo n-borne Infra 

theoretical models for PDRs fTielens fc HollenbachHl985j Red Telescone THIRTY JShihal etltl Il99lh and Balloou- 



ISternberg fc Dalgarncj [1989J) can explain the observed in- borne i n f rare d Carbon Explorer CBTCFi ijNakagawa, et all 



Il998j) have shown that the diffuse [CII] emission from the 
Send offprint requests to: B. Mookerjea lower density regions contributes significantly to the cool- 
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ing of the gas in our Galaxy and there is a significant 
drop in [CII] to continuum ratio towards regions of active 
star formation. Thus large scale high angular resolution 
studies of extended [CII] and CO emission from both the 
active star forming regions and PDRs with low densities 
and low UV radiation intensity are needed to understand 
the dominant cooling mechanisms under different physi- 
cal conditions within the molecular clouds. However the 
[CII] fine structure line is observable only from above the 
earth's atmosphere using aircraft and balloons, thus re- 
stricting the scope of the observational studies. Here we 
present the first results of a balloon-borne programme for 
observation of the [CII] emission from Galactic star form- 
ing regions. 

Situated at a distance of 450 pc, the Orion molec- 
ular cloud is the nearest region of active star forma- 
tion. It also show s clear evidence of ubiquity of PDRs 
1 Genzel fc Stutzkl Il989() . Owing to its proximity and 
brightness, numerous observational studies have been di- 
rected towards Orion and these have resulted in improving 
our understanding of the interplay between stars and the 
interstellar matter (ISM). In addition, Orion is also a good 
benchmark source to test any new observing system. The 
Orion molecular cloud also has the advantage of being 
like a ridge extended along the north-south direction with 
much smaller extent along the east- west. This is particu- 
larly useful for chopped observational studies, so that the 
emission lost due to self-chopping is negligible. 

Previous [CII] observations of Orion were carried 
out from the Kuiper Airborne Observatory. These in- 
clu de the mapp i ng of a 6.'5 x 10' region in OMC 1 
bv IStacev et alJ (|l993h . the mapping of 7'xl8' region 
in the Orion Molecu lar Ridge (OMC 1 and 2) by 
iHerrmarm et al.1 lll99jt) and heterodyne spec t roscop y of 
M42 bv lBoreiko et all frl988|k iBoreiko fc Betd l)l996l) . 

Here we present a large scale study of [CII] 158 /im 
fine-structure line emission towards Orion A. The map 
is fully sampled and covers an area of 30'xl5'. The rest 
of the paper is organized as follows: A brief description 
of the instrument, observational procedure and parame- 
ters are presented in Sect. 2. Sect. 3 presents the data 
analysis methods, Sect. 4 presents the results of observa- 
tion, Sect. 5 presents comparison with other observations 
and Sect. 6 presents a discussion of the intensity profiles 
and interpretation of the line intensity ratios using a PDR 
model. 

2. Instrument & Observation 

2.1. F PS 100 

A Fabry Perot Spectrometer (FPS) tuned to the [CII] 
157.7409 /jm line was installed at the foca l plane of the 
100 c m TIFR balloon-borne FIR telescope (jGhosh et all 
Il988|) . We refer to this new combination of instruments 
as FPS100. The FIR telescope has an f/8 Cassegrain 
configuration in which the secondary mirror can be wob- 
bled with an amplitude of 4f5 and a frequency of 10 Hz. 



The FPS was designed at ISAS, Japan and in an earlier 
configuration wa s used at the focal plane of the BICE 
l(Nakagawa et all Il998h . Interfacing of the FPS to the 
TIFR telescope necessitated modification of the optical 
components in order to achieve f-number compatibility be- 
tween the two. New electronics was also designed and built 
to implement the observing modes (discussed in Sect. l2.2|) . 
We point out that the FPS on the BICE had an angular 
resolution of 15', while on the TIFR telescope the achieved 
angular resolution is 1.'5. The FPS, in brief is a tandem 
Fabry Perot Spectrometer, consisting of two Fabry Perot 
interferometers: one a high order interferometer with mov- 
able plates to scan the wavelength (Scanning Fabry Perot; 
SFP) and the other a low-order interferometer with fixed 
plates (Fixed Fabry Perot; FFP) which acts as an order 
sorter for the SFP. The two interferometers together with 
other optics of the spectrometer and the FIR detector, are 
cooled to 2 K with liquid helium. T he FIR detector used is 
a stressed Ge:Ga photoconductor l|Hiromoto et all Il989() 
with an NEP of 3.1 x 10~ 15 Watt/VSz. 

2.2. Observation Modes 

The FPS100 has been designed for two observation 
modes, viz., chopped and unchopped. Both modes involve 
frequency modulation of the signal achieved by sweeping 
the SFP up and down at a preselected scan frequency and 
over a preselected velocity range around the expected po- 
sition of the [CII] line. The spectral scanning, in addition 
to being important for background reduction is also use- 
ful for estimating the neighboring continuum reliably. The 
primary difference between the two modes of observation 
lies in the fact that, in the chopped mode the sky chopping 
(effected by wobbling the secondary mirror) is functional, 
while in the unchopped mode there is no sky chopping. In 
the chopped mode the frequency of chopping is 10 Hz. For 
both modes it is possible to carry out spatial raster scans 
covering a rectangular box around the target source. 

Thus, the chopped mode involves sky chopping, spa- 
tial and spectral scanning. For meaningful measurements 
the frequencies of spatial scan (f sp ), sky chopping (vch) 
and spectral scan (v wv ) in this mode should maintain the 
following relation: u sp <^ W v<^ch- The signal detected in 
this mode is the difference between the emission in the 
two beams and is processed through linear phase band 
pass filters and further amplified (with selectable gains). 
There after the signal is phase sensitively detected (with 
online phase adjustment capability) and digitized for the 
telemetry. One of the outputs of the PreAmplifier for the 
SPcctrometer (PASP) is used as the output signal for the 
chopped mode. The signal is sampled at a frequency of 
10 Hz. 

In the unchopped mode, since there is no onboard 
background subtraction, the integration time over each 
spectral element is kept small, i.e., in effect a series of 
narrow band observations are done. Thus for this mode 
the spectral scan frequency is at least 4 times higher than 
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that for the chopped mode. For this mode the relation be- 
tween the different frequencies involved is z^ c h <v sp <i' wv . 
Since there is no modulation of the signal due to sky chop- 
ping (i/ c h = 0), the corresponding output from the PASP 
is directly digitized using a 16 bit ADC. The frequency of 
spectral modulation being higher in this mode, the signal 
is sampled at a frequency of 40 Hz. 

2.3. Observation Details 

The FPSIOO payload was flown for the first time from 
the TIFR Balloon Facility Hyderabad, in Central India 
(Latitude = N17?47, Longitude = E78?57) on 1999 
November 25. The observations presented here were car- 
ried out in the chopped mode. The choice of the chopped 
mode over the unchopped mode was based on considera- 
tions related to the noise level during the flights. For the 
chopped mode the noise level was found to be consistent 
with the pre-flight tests in which the flight condition was 
simulated by closing the window of the FPS and keeping 
the rest of the configuration identical to that during the 
flight. 

Jupiter was observed for absolute flux calibration as 
well as for determining the point spread function (PSF). 
The angular resolution achieved during these observations 
was 1 '5. A rectangular region of 30' x 15' towards Orion A 
was observed in approximately 50 minutes. The spatial 
scan rate was O.'32/s and the spectral scan frequency was 
0.5 Hz. 

3. Data Analysis 

The entire dataset for Orion A was split into two parts 
depending on whether the SFP is scanned from lower to 
higher frequency (Up) or in the reverse sense (Down) . The 
data analysis procedure described in the following text was 
applied independently to these two datasets. 

The signal measured is a combination of the emission 
from the astronomical source and the residual (post sky 
chopping) foreground radiation contributed by the atmo- 
sphere and the instrument. In the present work we have 
estimated the foreground emission by considering the sig- 
nal at the edges of the spatial scan lines where there is no 
apparent contribution from the astronomical source. As 
expected the residual emission in the chopped mode was 
found to be negligible. 

Each spectral scan was corrected for foreground radi- 
ation. Since the instrumental profile is much wider than 
the intrinsic astronomical line profile, the observed spec- 
tral profiles are dominated by the instrumental profile. 
Since the Lorentzian profile is a good approximation for 
the instrumental profile, we have fitted each spectral pro- 
file by a combination of the Lorentzian profile and a linear 
function of the wavenumber ct, representing the contribu- 
tion of the continuum emission at these wavelengths. The 
fitted function has the following form: 

f=~, !VI0V +C 2 o- + C 3 (1) 

2tt (cr - fi ) 2 + (T/2y 



The fitting procedure involves 5 parameters viz., Ci, 
C2, C3, r and piQ. piQ and T are respectively the central 
position and FWHM of the line. The observed line pro- 
file is slightly modified from the instrumental profile be- 
cause the astronomical line emission slightly broadens its 
width and shifts its central position as well. Thus for re- 
gions with strong [CII] line emission, all the 5 parameters 
are kept free and and T are determined accurately. 
Since th e [CHI linewid t h in t he region near 6 1 Ori C is 
5 km/s l|Boreiko et all Il988j) . Hence the line being not 
resolved by the FPSIOO we have kept fio and T fixed 
for all spectral profiles and determined the values of the 
remaining parameters from the fit. Fig.^presents the ob- 
served data (shown as points) corresponding to a spectral 
profile and the function (shown as a continuous line with 
Ato=63.404 cm" 1 and T= 0.045 cm" 1 ) fitted to it. This 
value of T corresponds to a velocity resolution of 210 km/s. 




1 1 1 1 1 

63.35 63.4 63.45 

Wavenumber (cm -1 ) 



Fig. 1. Typical Spectral Profile from Orion Observations. 
The observed data are shown as points and the function 
(combination of a Lorentzian and a straight line) fitted 
to the observations is shown as a continuous line. The 
fitted function corresponds to /io = 63.404 cm -1 and F = 
0.045 cm- 1 . 



This profile fitting results in two datasets correspond- 
ing to the intensities of the [CII] line and the continuum. 
Using the telescope aspect information the intensity (line 
or continuum) is gridded into a two dimensional matrix 
with pixels of size 0'A x 0.'4. The sky matrix thus derived 
is a convolution of the intensity distribution of the source 
and the bipolar (due to sky chopping) PSF. Using the PSF 
determined from Jupiter observations the data was decon- 
volved using a Maxim um Entropy Method (MEM) based 
deconvolution scheme l)Ghosh et all Il988|) to obtain the 
intensity distributions of the [CII] line and the continuum. 
The deconvolved maps of the [CII] line (and continuum) 
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intensity distributions obtained from the two independent 
datasets (Up and Down) match well. This confirms the 
reliability of most of the observed features. The final [CII] 
line and continuum maps are obtained by averaging the 
respective Up and Down datasets. 



4. Observational Results 

4.1. The [CII] map 

Fig.|21shows the fully sampled maps of velocity-integrated 
[CII] intensity and the continuum emission at 158 /im. The 
[CII] emission peaks at the position of 1 Ori C. Features 
detected in the [CII] emissi on are similar to those ob - 
served bv lStacev et all (1199.-^1 and lHerrmann et al 1 lll997l) . 
Both of these previous observations done with an effec- 
tive beamsize of 55" were confined (particularly along the 
right ascension) to the immediate vicini t y of t he central 
part of Orion A, while iHerrmann et alJ (l997|) extended 
the map northward along the ridge. The present observa- 
tions, though with slightly lower angular resolution cover 
a larger extent and detect a strong emission far towards 
the east as well as from M43 in the north. The [CII] map 
also exhibits an additional emission peak to the west of 
the main structure in the Orion A region. The extended 
nature of the [CII] emission supports the present under- 
standing that the [CII] emission not only arises from the 
direct interface between the ionized region (produced by 
the 9 1 Ori C) and the molecular cloud, but also from neu- 
tral atomic interfaces deeper into the molecular cloud pro- 
duced by the diffuse UV radiation field. There is no en- 
hancement of [CII] emission close to the outflow associated 
with BN/KL; this is consi s tent w ith the theoretical mod- 
els llHollenbach fc McKeeL Il989l) which predict that not 
much [CII] is produced in shocks. The [CII] emission falls 
off sharply towards the west, but appears to be more ex- 
tended in the east. This conforms w ell with the currentl y 
understood structure of the region (IWilson et alJ.ll997l) . 
This model assumes that the dense gas near the ionizing 
star shields the more distant gas to the east from the flux 
of the ionizing photons and to the west there being less 
dense neutral matter, the ionized gas flows in and domi- 
nates. The separated emission peak towards the west has 
also been detected in [CI] (1-0) and CO(3-2) emissions, as 
will be discussed in detail later. [CII] emission detected 
from the HII region M43, possibly arises from the associ- 
ated HII region which is energized by the star NU Ori. 



4.2. Comparison of [CII] & Continuum emission at 

158 \im 

Comparison of the [CII] emission with the continuum 
emission at 158 (right panel of Fig. |2J reveals that 
although there is an overall correlation between the two, 
there are distinct differences as well. The salient difference 
in features between the continuum and the [CII] emission 
are: (1) the positions of the peaks of the two emissions are 



different and (2) the [CII] emission is more extended than 
the continuum emission. 

In contrast to the [CII] emission which peaks at the 
position of 6 1 Ori C, the continuum emission peaks at the 
position of the BN/KL object. This possibly indicates that 
while the BN/KL object corresponds to the strongest em- 
bedded energizing source, giving rise to the dust emission 
peak, it is still not evolved enough to photodissociate gas 
which can emit [CII] strongly. An alternative or additional 
explanation for the difference in the peak positions could 
be the higher dust density close to the BN/KL object. 

The [CII] emission is more extended compared to the 
continuum emission, indicating the extended presence of 
neutral photodissociated gas and leaking of the UV pho- 
tons from 9 1 Ori C. The continuum emission in contrast 
falls off sharply with increasing distance from the embed- 
ded young stellar objects responsible for heating the dust 
in this region. 

However we point out that for the lower contour levels, 
the sensitivity of the continuum emission presented here 
is limited. This also accounts for the lack of the northward 
and eastward extension of the continuum emission as com- 
pared to the [CII] emission. This becomes apparent if we 
compare the features of the continuum map with those 
in the 138 fj,m map bv lMookeriea et all <|2000h . While the 
central features match exceedingly well, the extended na- 
ture of the continuum emission is clearly missed by the 
present observations. 



4.3. Gas Heating Efficiency 

The total [CII] luminosity detected in the central re- 
gion enclosed within contours up to 20% of the peak 
line emission is 45 Lq. The total FIR continuum lumi- 
nosity was calculated from the flux detected at 158 /im 
up to 6% of the peak, us ing the formula proposed by 
iThronson fc Telescol (|l986j) The calculation assu mes a 
dust temperature of 50 K (jMookeriea et al 1 hood and a 
dust cmissivity of A -1 and the calculated FIR continuum 
luminosity is 1.1 x 10 5 Lq. The choice of the lowest contour 
levels in the [CII] and continuum maps, up to which we 
integrate ensures integration over similar regions. Thus 
the [CII] luminosity detected is 0.04% of the total FIR 
continuum luminosity. 

Large scale survey of the Galactic [CII] emission 
()Nakagawa et all Il998l) using the BICE shows that the 
[CII] /FIR ratio is low towards the major star forming re- 
gions and the rati o is ~ 0.2% for c ompa ct sources. The 
ratio presented in iNakaeawa et alJ l)l998|) uses the total 
continuum emission derived between 40 and 120 //m. If 
the same ratio is calculated over 1 to 500 ^m, for a dust 
temperature of 50K and an emissivity index of 2.0 the ra- 
tio would be 0.1%. The BICE survey was conducted with 
an effective beamsize of 15', which is much larger than the 
resolution of our observations (1.'5). The resulting differ- 
ence in the beam filling factors of the two observations 
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Continuum at 158 /im 
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Fig. 2. Intensity maps of: Left : [CII] at 158 /im and Right :dust continuum (overlayed with the [CII] map at 158 /im 
contours in dotted lines). Peak of the [CII] intensity map is 3.9 x 10~ 3 erg s" 1 cm" 2 sr" 1 and that of the continuum 
is 4.76 x 10" 3 erg s" 1 cm" 2 sr" 1 . The contour levels are 95%, 90% to 10% (in steps of 10%) 5% and 2.5% of the 
respective peaks. Marked with the filled triangle is BN/KL ( 05^ 32™ 46?7 -05° 24' 17"(1950)), * is 9 1 Ori C ( 05^ 
32™ 49?0 -05° 25' 16"(1950)) and the filled square is M 43 ( 05 /l 33™ 03?0 -05° 18' 06"(1950)). 



can possibly account for this discrepancy in the estimated 
ratio of the [CII] and the FIR luminosities. 

The reduction in the [CII] to continuum ratio towards 
regions of high mass star fo rmation is explained by model 
calculations. PDR models jHollenbach et all Il99ll) show 
that in high density regions with large Go, the [CII] line 
is rather saturated due to its lower excitation energy and 
lower critical density, while lines such as the [O I] 63 /mi 
line become dominant coolants. In addition, dust grains 
being positively charged in regions with high Go, the effi- 
ciency of photoelectric heating of gas decreases. 



5. Morphological comparison with other 
observations: 

We present a comparison of the observed [CII] and dust 
continuum emission with the 13 CO(1-0), [CI] (1-0), CO (3- 
2) and the radio continuum emission at 1.5 GHz. The 
comparison with 13 CO(1-0) is aimed at understanding the 
overall column density structure of this region, vis-a-vis 
the presence of heating sources. The comparison with [CI] 
and CO (3-2) emissions is aimed towards understanding 
the volume density structure and the distribution of the 
UV field, in an attempt to derive a coherent picture of 
PDRs giving rise to all of these emissions. The compari- 
son of [CII] emission with the radio continuum emission is 
intended to identify ionized gas contributing substantially 
to the [CII] emission. 




5 h 33 ra 20 5 4 0*0 

RA (1950) 



Fig. 3. Overlay of the [CII] at 158 /im (con tours as dotted 
lines) and 13 CO (1-0) jPlurne et all 12000^ emissions from 
Orion A. 13 CO (1-0) has been smoothed to a resolution of 
1.'5. Peak 13 CO (1-0) intensity is 32.7 K km s _1 , contours 
levels are 10% to 90% in steps of 10% and 95% of the 
peak. The contour levels for the [CII] map are 2.5 10 30 
50 70 90 % of the peak intensity. The marked sources in 
the map are identical to Fig. [5] 



5.1. Comparison with 13 CO(1-0): Column density 
tracer 

Observations sugge st that 12 CO is optically thick 
l|Castets et allll990|) in Orion A. On the other hand, 13 CO 
is optically thin in most cases and is a better tracer of the 
column density of the neutral molecular gas. We compare 
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our [CII] data with that of 1 3 CO(1-0) observed with a 
beam of 47", at the FCRAO bv lPlume et all (l200Ctl . Fig.0 
also shows an overlay of the [CII] emission with the 13 CO 
(1-0) emission. The 13 CO (1-0) data has been integrated 
over velocities ranging between 5 and 12 km/s. smoothed 
to the resolution of the [CII] map (1.'5) and regridded on 
a 0.'4 grid. 

The velocity integrated 13 CO (1-0) intensity map 
shows two prominent peaks, one coinciding with the in- 
frared object BN/KL and the other lying to the south of 
it. The continuum emission also peaks at the location of 
the BN/KL object, thus implying a core of high column 
density (and possibly high volume density as well) with 
an embedded heating object at a reasonably early stage of 
evolution. The southward extension of the dust continuum 
emission is resolved into a clump in the 13 CO(1-0) map. 
However, neither the 13 CO(1-0) nor the dust continuum 
map show the separated emission peak to the west, which 
is seen in the [CII] map. The [CII] emission in contrast 
shows a single peak at the position of the ionizing star 9 1 
Ori C. Since both [CII] and 13 CO (1-0) emission are to a 
good extent optically thin, the dissimilarity in the overall 
emission features hint at the difference in the phases of the 
ISM which dominantly contribute to these emissions. The 
13 CO(1-0) emission originates from cold molecular mate- 
rial and hence is more widespread. We note that there 
is a distinct similarity in the eastward extension of the 
13 CO(1-0) emission with that of the [CII] emission. The 
13 CO(1-0) map shows an isolated peak at the position of 
M43 to the north, which is a little offset from the [CII] 
intensity enhancement in this direction. 

5.2. Comparison with Cl(l-O), CO(3-2): PDR view 

The present understanding of the PDRs in molecular 
clouds suggests that the [CII] emission from the PDRs 
arises from the surface of UV irradiated clumpy molec- 
ular clouds. The penetration of UV into the molecular 
material creates a stratification in which several carbon 
bearing species appear and disappear sequentially with 
increasing visual extinction into the cloud. Other species 
which are present in the PDRs and contribute to the cool- 
ing of these regions include, the neutral atomic carbon 
and the CO (mid-J rotational transitions) from close to 
the edge of the UV irradiated surface and the low- J rota- 
tional transitions of CO and its isotopomers which trace 
the large scale gas distribution of the molecular cloud. 

Fig.Hshows the measured [CI] 3 Pi- 3 P {left) and CO 
(3-2) {right) emission from the same region of Orion A, 
as we have mapped in [CII]. The [CI ] and CO (3-2) dat a 
have been taken from the paper bv llkeda et al.l l|l999h . 
The [CI] and the CO (3-2) data were observed with an- 
gular resolutions of 2'.2 and 3' respectively. Both datasets 
were presented in the paper with an angular resolution of 
3' and on a U5 grid spacing. For proper morphological 
comparison, we have smoothed the [CII] data to a resolu- 
tion of 3'. 



[CI] arises only from the neutral phase of the ISM, 
while CO (3-2) traces the warm and/or dense gas. In 
other words the CO (3-2) emission stems from the pos- 
sible clumpy structures in the ISM which also harbour 
the [CII] and [CI] emitting regions closer to their sur- 
faces. The [CI] emission originating from the neutral gas 
is more widespread than the [CII] emission and peaks at 
two positions, none coinciding either with the BN/KL ob- 
ject or 6 1 Ori C. The southern peak matches with the 
second peak of the 13 CO(1-0) map. We point out that 
the peaks of the emissions of the three lines follow the 
sequence [CII] /CO/ [CI], which is different from the PDR 
scenario in which the dominant carbon bearing species 
changes from C + to C° to CO. The eastward extension 
and the separated emission peak to the west seen in the 
[CII] map are clearly seen in the [C I] map also. We note 
that [CI] emission shows no enhancement at the position 
of M43 contrary to the expectation that presence of [CII] 
emission indicates the presence of PDR and hence [CI] 
emission. This indicates that the [CII] emission in M43 
originates from the associated HII region, rather than from 
the neutral PDR. 

At this angular resolution the CO (3-2) map shows 
much less features as compared to the [CI] map. The emis- 
sion peaks at a single position centered on the BN/KL ob- 
ject, much like the continuum emission. The peak position 
is identical to the northern peak of the 13 CO(1-0) emis- 
sion. The extension towards the east and the separated 
emission peak to the west as seen in the [CII] map arc also 
clearly distinguishable in the CO(3-2) map. Presence of 
these features in the CO(3-2) map suggests enhancement 
of temperature and hence presence of embedded sources. 
This is also substantiated by the observed [CII] emission. 

5.3. Comparison with 1.5 GHz radio continuum 

The first ionization potential of carbon is 11.26 eV (less 
than that of hydrogen). Hence, the [CII] emission stems 
not only from the neutral PDRs but also from the ion- 
ized gas in the HII regions. Model calculations of dust 
free HII regions imply that stars hotter than spectral type 

07 produce sufficiently high flux of photons with energies 
greater than 24.587 eV to maintain all of the helium in 
the form of He + , and therefore all of the carbon in the 
form of C ++ (the ionization energy required to produce 
C ++ lies only 0.204 eV below that required to ionize He). 
However, presence of dust within the HII region causes 
the stellar Ly continuum radiation field to become de- 
pleted of He-ionizing photons. As a result some C + may 
be found in dusty HII regions ionized by stars hotter than 
07. Alternatively, the HII regions whose ionizing stars are 

08 or cooler are incapable of maintaining a significant 
fraction of carbon in the C ++ state and hence may have 
a significant population of C + . 

Fig-Elshows an overlay of the 1.5 GHz radio continuum 
emissio n measured using VLA , with an angular resolution 
of V bv ISubrahmanvan et al ] l|200ll) . with the [CII] map. 
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Fig. 4. Overlay of [CII] map (with Left : [CII (1-0) fc Right : CO (3-2) maps of Orion A (Resolution of 2.'1 & 3' 
respectively, CI data smoothed to 3': llkeda et ailll99flh . Peaks of [CI] & CO(3-2) arc 41.4 K km s" 1 and 382.4 K km 
s _1 respectively. Contour levels for both maps are 10% to 90% in steps of 10% and 95% of the peak. Sources marked 
are the same as in Fig. [2 Both maps are at an angular resolution of 3' and gridded on grid. Relevant details for 
the [CII] map as well as the marked sources in the map are identical to Fig. |3J 



-5°12'00" - 
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-5°28'00" - 




-5°36'00" - 

i I i I i L_ 

5*33"Z0?Q 40"0 00"0 

RA (1950) 

Fig. 5. O verlay of VLA 1.5 GHz co ntinuum image of 
Orion A ijSubrahmanvan et"aT1. l200l|) with the [CII] in- 
tensity map. The peak intensity is 26 Jy beam -1 and the 
contours are at 0.25, 1.0, 2.5, 5, 10, 30, 50, 70 and 90% of 
the peak. Relevant details for the [CII] map as well as the 
marked sources in the map are identical to Fig. 

The basic features of the radio continuum emission arc 
very different from the [CII] emission, implying origin in 
different components (ionized and/or neutral) of the ISM. 
The [CII] emission is extended along the north-south di- 
rection, while the radio continuum emission is extended in 
the east-west direction. The [CII] halo is more extended 
toward the east, while the radio continuum is extended 
mostly toward the west. This agrees well with the recent 



models for the Orion A region {Wilson et all Il997|) . The 
radio continuum emission peaks at the positions of the two 
HII regions associated with M42 and M43. The [CII] map 
peaks at the position of M42 and shows enhanced emis- 
sion close to M43, with no corresponding enhancement in 
[CI] emission. 

We note, M42 is ionized by a 06-07 star, so that the 
calculate d [CII] intensity orig inating in the ionized gas is 
only 5% {Russell et all ll980|l . However M43 is ionized by 
a BO. 5 V star, so the enhancement of [CII] emission with 
no corresponding [CI] emission peak can be attributed to 
the ionized gas in the HII region. 

6. Discussion 

6.1. Line Intensity Profile &. Ratios 

In order to probe the density and the UV field strength 
distribution with increasing distance from 6 1 Ori C, we 
derive profiles of the [CII], [CI], CO (3-2) line intensities 
for a east- west cut at the declination of the 1 Ori C and 
calculate the intensity ratios at a few selected points. All 
three datasets have/are smoothed to an angular resolution 
of 3' and are gridded on a 1.'5 grid. 

Fig.Elshows the intensity profiles referred to above, of 
[CII], [CI] and CO (3-2) along the east- west direction at 
the declination of 9 1 Ori C. The [CII] emission falls off 
rapidly (steeper towards the west) away from 9 1 Ori C, 
but continues showing residual emission well beyond. The 
[CI] emission arising from cold neutral medium is more 
widespread. The CO (3-2) emission tracing warmer gas is 
moderately distributed, and it's profile shows a hump to- 
wards the east which is similar to that seen in the [CII] 
intensity profile. 
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Fig. 6. East-west profile of intensities of [CII], [CI] and CO(3-2) at the declination of 1 Ori C. The [CI] and CO(3-2) 
intensities have been multiplied by a factor of 100. 



Table ^ presents the absolute [CII] intensities, the C + 
column densities and ratios observed at several selected 
positions (offsets relative to 6 1 Ori C) along the E-W cut 
shown in Fig.EJas well as the position of the BN/KL object 
(-90", 90"). In order to obtain a first order estimate the 
C + column densities have been calculated assuming the 
line to be optically thin, the densities to be high enough 
so that the emission is th ermalized to a temperature of 
165 K ijStacev et all Il993f) and the beam filling factor to 
be unity. The C + column densities are thus lower limits. In 
Tabled the UV flux and density correspond to the results 
of the PDR model analysis of these intensities. We discuss 
the PDR modeling in details in the next subsection. 



6.2. Density and UV Flux from PDR Model 



We have used the PDR model of lKaufman et all l|l999j) 
to interpret the observed intensity ratios of [CII], [CI] 
and CO (3-2) lines and to estimate the volume density 
and the UV field strength at selected positions. Tabled 
presents the observed line ratios and results of model- 
ing including the values of t he reduced \ 2 at eac h of 
these positions. The models bv lKaufman et all l)l999(l are 
similar in terms of chemistry, radi ative transfer, thermal 
balance etc. to those presented bv iTielens fc Hollenbachl 
l)l98^ . lHo11enba,cb et a,l.Ml99lh and IWolfire et a,l]|jl99(t . 
But they use updated reaction rates, heating rates as well 
as chemical abundances and most importantly they take 
into account the heating due to PAHs which is important 
for the i nterpretation of the F IR cooling lines like [CII] 
and [OIJ. lKaufman et al.ll|l999l) have created a database of 



PDR models, for a density range from n= 10 to 10 cm 
and for a far-ultraviolet (FUV) radiation field ranging be- 
tween Go = 0.3 to 3 x 10 6 (where Go is in units of an 
average interstellar flux of 1.6 x 10 -3 ergs s _1 cm -2 ). 

In order to determine the densities and UV fluxes 
which can explain intensity ratios, we have done a x 



fitting of the observed ratios with respect to the model 
predictions. Fig. shows in grey scales the reduced \ 2 
for the observed line ratios at the selected positions. The 
grey scale is proportional to the \ 2 an d weighted with the 
individual observational uncertainties. Based on the esti- 
mated calibration errors of the individual observations, we 
take an uncertainty of 30% for the ratios involving [CII], 
while for the [CI]/CO(3-2) ratio the uncertainty is taken 
to be 20%. The lines correspond to the observed intensity 
ratios and the absolute [CII] intensity (as described in the 
legends). The values for the volume density and UV flux 
at each position were determined from the location of the 
minimum \ 2 ■ We discuss the result of modeling for each 
position separately. 



6.2.1. Positions away from the [CII] peak 

(-90" ,90" ): This position corresponds to the peak of the 
CO (3-2) map and is approximately the location of the 
BN/KL object. The best fit model corresponds to an UV 
field of Go = lO 4 25 and a volume density of 10 5 cm -3 . 
The density compares well wit h those derived from the 
multi-line study of Orion A bv IWilson et all i|200ltl . We 
point out that the fitted model is not accurate; though 
it explains the line intensity ratios, the derived physical 
conditions predict an absolute [CII] intensity lower than 
the observed value. 

East of 9 1 Ori C: We are able to explain all the inten- 
sity ratios as well as the [CII] intensity at the 2 eastern 
positions (180",0") and (360",0") using the PDR models. 

West of 1 Ori C: We restrict the analysis of the west- 
ern side to a single point (-180", 0") owing to the limited 
dynamic range of the [CII] map. The model fitted is ac- 
ceptable, although we note that the predicted [CII] inten- 
sity is lower than the observed intensity. 



B. Mookerjea et al.: 158 /im [CII] emission from Orion A 9 



Table 1. Ratios of observed line intensities (in K km s 1 ) and the modeled parameters, at selected positions. 



Position 1 


\CII) 
ergs s _1 cm -2 sr _1 


N c+ 
cm -2 


[On] 
"[erf 


[Crl] 

CO(3-2 


^ [Cr] 

> CO(3-2) 


Density 
log(cm -3 ) 


UV flux 

log(Go) 


x 2 


Remarks 


(-90",90") 


1.30(-3) 


7.3 (17) 


267 


78 


0.29 


5.0 


4.25 


1.0 


CO(3-2) peak 


(0",0") 


2.55(-3) 


1.4 (18) 


631 


220 


0.35 


5.25 


5.5* 




[CII] peak 


(180",0") 


8.57(-4) 


4.8 (17) 


361 


141 


0.39 


5.0 


4.75 


0.7 




(-180",0") 


4.21(-4) 


2.4 (17) 


84 


30 


0.35 


4.75 


2.5 


1.6 




(360",0") 


5.30(-4) 


3.0 (17) 


239 


104 


0.43 


4.75 


3.75 


0.6 




1 Offsets relative to the position of 


9 1 Ori C (a 


1950 = 


05 h 32 m 


49:0, S 1950 = 


-5°25'16") 









* Since the PDR model gives no solution, the UV flux has been calculated using the lu minosity of Trap ezium cluster (L~ 10 Lq) 
and the distance between the Trapezium cluster and the molecular cloud (d = 0.25pc lO'Dell I \2QQ ll) ) 



6.2.2. The [CII] peak 

((/',(/'): This corresponds to the position of the star 9 1 
Ori C. For the parameter space explored the models do 
not have any solution. We identify the failure as due to the 
excess of observed [CII] emission over the model predic- 
tions. We note that if models corresponding to higher UV 
field were considered, a solution could be derived, but it 
would not have been realistic. This is because for the given 
luminosity (~ 10 5 Lq) of the ionizing stars, the Trapezium 
cluster , and their dis tance of 0.25 pc from the molecular 
cloud 'Dell I l200l[) the strength of the UV field at this 
position would at most be 3.3 x 10 5 Go- 

Since the ionizing star is of the type 06 the contri- 
bution of the ionized gas towards the [CII] emission will 
be negligible, unless the high dust density plays a crucial 
role. However this scenario for the present is highly spec- 
ulative, so we do not persist with it. Instead we keep in 
mind that the geometry and structure of Orion A is far 
too complicated and not fully understood to be satisfacto- 
rily modeled by the presently available PDR models. We 
consider the effect of inappropriate model geometry as a 
probable reason behind the failure of the PDR models to 
generate a satisfactory solution. 

We first estimate the C + column density which is 
needed to explain the detected [CII] intensity, assuming 
that the line is optically thin. For regions outside HII re- 
gions the almost temperature independent, critical density 
(n cr j t ) for collisional excitation is 4 x 10 3 cm -3 . Since the 
measured volume densities for these regions are well in ex- 
cess of these values (typically ~ 10 5 cm~ 3 ) the [C II] emis- 
sion is in thermal equilibrium. IStacev et alJ l|l993|) suggest 
that the minimum gas kinetic temperature of the region is 
~ 165 K. With these assumptions and values we estimate 
the C + column density to be 1.4 x 10 18 cm~ 2 . Assuming 
an abundance ratio of [C + ]/[H] ~ 3 x 1CP 4 , this corre- 
sponds to a hydrogen column density of 4.8 x 10 21 cm~ 2 , 
which is equiva l ent to an A v of 2.5 mag. The model of 
iKaufman etaTI ljl99fl|) treats the PDRs for A v upto 10 
mag, but being a homogeneous slab model considers the 
[CII] emitting region to be confined to the surface. If in- 
stead of a single such surface, a number of such surfaces 
(optically thin for [CII] emission) are placed, then it is 
possible to enhance the C + intensity relative to the other 
line intensities. The need to invoke a number of such sur- 



faces, suggests the presence of non-homogeneity or dumpi- 
ness in the ISM. The scenario changes appreciably if a 
spherical clump model of PDRs is used, in which case the 
width of the [CII] e mitting layer chang es depending on the 
size of the clump (|Koster et allll9*9*^l . A model treating 
an ensemble of spheres is more realistic, since observa- 
tions support that the ISM consists of clumps immersed 
in a tenuous inte r-clump medium, rathe r than being a 
homogeneous slab. lMochizuki et alJ l)l994j) have explicitly 
discussed that in the clumped scenario, the [CII] emis- 
sion is enhanced significantly compared to the homoge- 
neous case, if the average size of the clumps are such that 
A v < 10 mag. 

6.3. Summary of PDR modeling 

The PDR models are able to explain the observed line 
intensity ratios, satisfactorily for positions away from the 
energizing sources like the 8 1 Ori C and the BN/KL ob- 
ject. At the position of the [CII] peak in particular it fails 
to provide a solution, mainly because of the excess in [CII] 
emission over the model prediction. This is probably be- 
cause the actual geometry is far more complex than the 
model assumptions. However, in the absence of velocity- 
resolved [CII] observations it is difficult to throw more 
light on the actual structure of the [CII] emitting region. 
We note that although the contribution of the ionized gas 
towards the [CII] emission at the peak position in M42 
is expected to be negligible it may not be totally absent 
and must be taken into account while attempting better 
constrained modeling at other positions. 

According to the PDR model, the [CI] /CO (3-2) in- 
tensity ratio constrains the volume density, while the 
[CII] /[CI] is very sensitive to the strength of the UV flux. 
For all the selected positions, the [CI]/CO(3-2) intensity 
ratio varies very little and the derived volume densities 
are almost the same. The larger variation of [CII] /[CI] in- 
tensity ratios corresponds to a larger spatial variation of 
the UV field strength. 

7. Summary 

We have presented the first results of observations with a 
Fabry Perot spectrometer tuned to [CII] 157.7409 ^m on- 




log 10 n [cm 3 ] 

Fig. 7. Each panel corresponds to a selected position, the offsets of which are expressed relative to the position of 6 1 
Ori C. Plotted in grey scale is the reduced-% 2 (absolute values shown in the wedge) calculated from a comparison of 
the predicted intensity ratios and [CII] intensities by the Kaufman model, over the entire parameter space of volume 
density (in units of cm~ 3 and UV flux (in units of Habing field). The darker regions correspond to lower values of the 
reduced-x 2 - The details of the drawn lines are given in the top left panel and are identical for all panels. 



board a balloon-borne FIR telescope, FPS100. These ob- 
servations demonstrate the capabilities of the instrument 
we had planned for. The line integrated intensity map of 
[CII] emission from the Orion A region as well as the con- 
tinuum emission at 158 /xm have been obtained here. The 
[CII] luminosity is 0.04 % of the FIR continuum luminos- 
ity. Morphological comparison of the [CII] map with with 
[CI] at 492 GHz, CO (3-2) and 13 CO (1-0) data show (i) 
that the peaks of the emission occur in the sequence of 
[CII] /CO/ [CI] which is contrary to the PDR scenario and 
(ii) the [CII] emission from M43 in the absence of corre- 
sponding [CI] emission stems from the HII region and not 
from the PDR. 

The [CII], [CI] and CO (3-2) line inte nsity ratios have 
been modeled using the PDR models of iKaufman et all 
l)l999l) . The models though somewhat simplistic, consid- 



ering homogeneous plane parallel slabs instead of the 
clumpy structure seen observationally, are a first step to- 
wards understanding the observations. The models do not 
reproduce the observed [CII] intensity as well as the inten- 
sity ratios involving [CII] for positions which are close to 
the ionizing source (viz., 8 1 Ori C and possibly BN/KL). 
This possibly is an effect of the assumption of plane par- 
allel geometry instead of the more complicated structure 
of the region. The models produce acceptable solutions 
for regions with diffuse UV radiation. Similar to the [CII] 
intensity profile both the incident UV flux and the vol- 
ume density show sharper decline towards the west of 6 1 
Ori C, than towards the east. For all the selected positions 
the volume density derived using the PDR models varies 
little, while the UV field strength varies substantially. 
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